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Abstract. We present the results of dust continuum and molecular line observations of two massive cluster-forming 
clumps, NGC 2264-C and NGC 2264-D, including extensive mapping performed with the MAMBO bolometer 
array and the HERA heterodyne array on the IRAM 30m telescope. Both NGC 2264 clumps are located in 
the Mon OBI giant molecular cloud complex, adjacent to one another. Twelve and fifteen compact millimeter 
continuum sources (i.e. MMSs) are identified in clumps C and D, respectively. These MMSs have larger sizes 
and masses than the millimeter continuum condensations detected in well-known nearby protoclusters such as 
p Ophiuchi. The MMSs of NGC 2264 are more similar in size to the DCO + 'cores' of p Oph, although they are 
somewhat denser and exhibit broader linewidths. Most of the MMSs of NGC 2264-C harbor candidate Class 
protostars associated with shocked molecular hydrogen jets. Evidence for widespread infall motions is found in, 
e.g., HCO + (3-2) or CS(3-2) in both NGC 2264-C and NGC 2264-D. A sharp velocity discontinuity ~ 2 km.s" 1 in 
amplitude is observed in N2H + (l-0) and H 13 CO + (1-0) in the central, innermost part of NGC 2264-C, which we 
interpret as the signature of a strong dynamical interaction between two MMSs and their possible merging with 
the central MMS C-MM3. Radiative transfer modelling supports the idea that NGC 2264-C is a highly unstable 
prolate clump in the process of collapsing along its long axis on a near free-fall dynamical timescale ~ 1.7xl0 5 
yr. Our model fit of this large-scale collapse suggests a maximum mass inflow rate ~ 3xl0 -3 Mq yr -1 toward 
the central protostellar object C-MM3. In NGC 2264-D, we estimate a mass infall rate Mdmmi ~ 1.1 x 10" 4 M 
yr -1 toward the rotating Class object D-MM1, also based on radiative transfer modelling of the observations. 
Such infall rates are sufficiently high to overcome radiation pressure and allow the formation of ~ 20 Mq stars by 
accretion in ~ 1.7x 10 yr, i.e., a time similar to the global dynamical timescale of the central part of NGC 2264-C. 
We conclude that we are likely witnessing the formation of a high- mass ( ^ 10 Mq) protostar in the central part 
of NGC 2264-C. Our results suggest a picture of massive star formation intermediate between the scenario of 
stellar mergers of Bonnell et al. (1998) and the massive turbulent core model of McKee & Tan (2003), whereby 
a turbulent, massive ultra-dense core is formed by the gravitational merger of two or more Class protostellar 
cores at the center of a collapsing protocluster. 

Key words, stars: formation - stars: circumstellar matter - stars: kinematics - ISM: clouds - ISM: structure - 
ISM: kinematics and dynamics - ISM: individual object: NGC 2264 



1. Introduction 

While most stars are believed to form in clusters (e.g. 
Adams & Myers 2001, Lada & Lada 2003), our present 
theoretical understanding of the star formation process is 
essentially limited to isolated dense cores and protostars 
(e.g. Shu et al. 1987, 2004). Detailed studies of the earliest 
phases of clustered star formation are crucially needed if 
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we are to explain the origin of the stellar initial mass func- 
tion (IMF) and the birth of massive stars. There is indeed 
a growing body of evidence that high-mass (M* > 8M ) 
stars may be able to form only in a clustered environment 
(e.g. Zinnecker et al. 1993, Testi et al. 1999, Clarke et al. 
2000). 

On the theoretical side, two main scenarios have been 
proposed for clustered star formation in turbulent molecu- 
lar clouds. In the first scenario, the distribution of stellar 
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masses is primarily determined by turbulent fragmenta- 
tion. Briefly, self-gravitating condensations (each contain- 
ing one local Jeans mass) form as turbulence-generated 
density fluctuations (e.g. Padoan & Nordlund 2002), then 
decouple from their turbulent environment (e.g. Myers 
1998), and eventually collapse with little interaction with 
their surroundings. In this scenario, high-mass stars are 
built-up by a scaled up version of the accretion process 
believed to be at work in low-mass protostars: a high 
accretion rate (> 10 -3 Mo.yr -1 ), generated by the high 
pressure of the turbulent environment (e.g. McKee & Tan 
2003) and/or the influence of an external trigger (e.g. 
Hcnnebelle et al. 2003), is sufficient to overcome the radi- 
ation pressure that would normally halt the collapse soon 
after M* w 8M Q (e.g. Stahler et al. 2000). 
By contrast, in the second scenario, the distribution of 
stellar masses results entirely from the dynamics of the 
parent protocluster. Here, individual protostellar seeds 
have large proper motions with respect to the mean cloud 
velocity. They accrete cloud material competitively while 
orbiting in the gravitational potential well of a collaps- 
ing protocluster and possibly collide with one another 
(e.g. Bonncll et al. 1998, 2001). In this alternative sce- 
nario, competitive accretion and dynamical interactions 
between individual cluster members play a dominant role 
(e.g. Klessen & Burkert 2000, Bate et al. 2003). In the 
dense central region of the collapsing protocluster, the 
probability of encounters becomes large enough that mas- 
sive stars can form by coalescence of intermediate-mass 
protostars (e.g. Bonnell & Bate 2002). Determining which 
of these two pictures dominates in actual cluster-forming 
clouds is a major open issue. 

As one expects the internal dynamics of protoclus- 
ters to differ markedly in the two pictures, comprehensive 
molecular line studies of the kinematics of prestellar con- 
densations and Class protostars in cluster-forming re- 
gions can provide powerful observational tests (cf. Andre 
2002). In an effort to improve our knowledge of the dy- 
namics of protoclusters and constrain theoretical models 
of clustered star formation, we carried out a combination 
of millimeter dust continuum and molecular line observa- 
tions of the active cluster-forming region NGC 2264, con- 
taining more than 360 near-IR sources (Lada et al. 1993, 
Lada & Lada 2003). Located in the Mon OBI molecu- 
lar cloud complex (d~ 800 pc), NGC 2264 is known for 
its famous Cone Nebula. Three arcminutes or ~ 0.7 pc 
to the north of the Cone Nebula, Allen (1972) discov- 
ered a bright embedded IR source, hereafter called IRS1, 
associated with IRAS 06384+0932, and also known as 
Allen's source. IRS1 is a B2-type object with estimated 
visual extinction Ay ~ 20-30 (Allen 1972, Thompson et 
al. 1998), bolometric luminosity Li, i ~ 2300 L© (Margulis 
et al 1989), and mass M* ~ 9.5 M Q (Thompson et 
al. 1998). Approximatively 6 arcminutes or ~ 1.4 pc to 
the north-west of IRS1, there is another IRAS source, 
IRAS 06382+0939, hereafter called IRS2. IRS2 is a Class I 
young stellar object (YSO) with L^i ~ 150 L© (Margulis 
et al. 1989). Both IRS1 and IRS2 are associated with 



molecular outflows, named NGC 2264-C and NGC 2264- 
D, respectively (Margulis et al. 1988), and dense molecu- 
lar clumps revealed by CS(2-1) (Wolf-Chase et al. 1995) 
and submillimeter dust continuum (Ward- Thompson et al. 
2000; Williams & Garland 2002, Wolf-Chase et al. 2003) 
observations. An H2O maser was found in each clump 
(Genzel et al. 1977; Mendoza et al. 1990), and a 44-GHz 
methanol maser was found in NGC 2264-C (Haschick et al. 
1990). Methanol masers are known to be clear signposts 
of ongoing intermediate/high-mass star formation (Minier 
et al. 2001). 

Previous kincmatical studies of NGC 2264-C have revealed 
a complex velocity structure (Kriigel et al. 1987). Evidence 
for large-scale collapse in the entire clump, disrupted by 
outflow motions on small scales, was reported by Williams 
& Garland (2002). Triggered star formation, due to an 
expanding shell ~ 0.12 pc in diameter around IRS1, was 
suggested by Nakano et al. (2003) based on observations 
with the Nobeyama millimeter interferometer. These pre- 
vious studies had either lower angular resolution and/or 
were focused on a much smaller region than the observa- 
tions reported here. 

The layout of the paper is as follows. Observational details, 
dust continuum results, and molecular line results are de- 
scribed in Sect. 2, Sect. 3, and Sect. 4, respectively. Section 
5 discusses the nature of the embedded millimeter sources 
identified in NGC 2264, while Sect. 6 and Sect. 7 present 
radiative transfer models of NGC 2264-C and NGC 2264- 
D-MM1, respectively. We compare our observational re- 
sults with various scenarios of high- mass star formation in 
Sect. 8. Our main conclusions are summarized in Sect. 9. 

2. Millimeter Observations 

We performed 1.2 mm dust continuum mapping obser- 
vations of NGC 2264 in December 2000 with the IRAM 
30 m telescope near Pico Veleta, Spain, using the MAMBO 
37-channel bolometer array. Eleven on-the-fly maps were 
taken in the dual-beam scanning mode, with individual 
sizes ranging from 13' x 10' to 5'x4'. Chopping was per- 
formed by wobbling the secondary mirror of the telescope 
at 2 Hz with a throw ranging from 32" to 72". The 
scanning speed was either 6"/sec or 8"/sec. The atmo- 
spheric opacity at zenith varied between 0.11 and 0.45 at 
A = 1.25 mm. Pointing and focus checks were made every 
hour. The pointing accuracy was measured to be better 
than 4" and the angular resolution was ~ 11" (HPBW). 
Calibration was achieved by mapping Mars (primary cali- 
brator) and is believed to be reliable to better than ~ 15%. 
Our resulting 1.2 mm continuum mosaic of the NGC 2264 
region covers a total area of ~ 400 arcmin 2 , i.e, ~ 23 
pc 2 with a mean rms noise a = 7 mJy/ll"-beam (i.e. 8 
mJy/13"-beam, see Fig.GJ. 

We performed follow-up molecular line observations 
with the IRAM 30 m telescope in March and May 
2002. The observed molecular transitions were N2H + (l-0), 
C 34 S(2-1), HCO+(1-0), H 13 CO+(1-0) at 3 mm, CS(3-2), 
C 34 S(3-2) at 2 mm, CS(5-4), C 34 S(5-4), HCO+(3-2) and 
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Fig. 1. Millimeter dust continuum mosaic of NGC 2264-C and NGC 2264-D smoothed to an effective angular resolution 
of 13". The first white contour level is at 30 mJy/13"-beam and corresponds to a column density ~ 2xl0 22 cm~ 2 
(assuming k=0.005 cm 2 .g _1 and Td—15 K). The other white contour levels go from 120 to 300 mJy/13"-beam by 90 
mJy/13"-beam, and from 300 to 1100 mJy/13"-beam by 200 mJy/13"-beam. The mean rms noise is ~ 8 mJy/13"- 
beam. The open white star symbols show the positions of the two IRAS sources IRS1 and IRS2. 



H 13 CO+(3-2) at 1.1 mm. The half-power beam width of 
the telescope was ~ 26", ~ 17", ~ 10"at 3 mm, 2 mm, 1.1 
mm, respectively. For most transitions, we took maps in 
the on-the-fly mode around the peaks seen in the 1.2 mm 
dust continuum map of Fig. ^ We used either four SIS het- 
erodyne receivers simultaneously or, for the 1.2 mm band, 
the HERA 9-pixel heterodyne array. As backend we used 
an autocorrelation spectrometer with a spectral resolution 
of 20 kHz at 3 mm, and 40 kHz at 2 mm and 1.2 mm. The 
corresponding velocity resolution ranged from 0.05 to 0.08 
km.s -1 per channel. 



3. Dust continuum mapping results 

3.1. Source extraction 

Our 1.2 mm dust continuum maps were reduced and 
combined with the IRAM software for bolometer-array 
data ("NIC"; cf. Broguiere et al. 1995). Combining all 
our maps we obtained the image shown in Fig. ^ The 
resulting mosaic confirms the presence of two distinct 
clumps, NGC 2264-C and NGC 2264-D, which closely 
follow the CS(2-1) integrated intensity map of Wolf-Chase 
ct al. (1995). 
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Table 1. Source extraction results in NGC 2264-C and NGC 2264-D 



Source" 


Coordinates 


Undec.FWHM" 


P.A. C 


(mJy/ll"-beam) 


S 


Slack 




Q2000 


$2000 


(arcsecond) 


(deg) 


(mJy) 


(mJy/ll"-beam) 


C-MM1 


06:41:18.0 


+09:29:02 


14x12 


16 


255 ± 11 


354 


186 


C-MM2 


06:41:15.3 


+09:29:10 


22x13 


-12 


183 ± 10 


433 


296 


C-MM3 


06:41:12.4 


+09:29:11 


18x13 


123 


573 ± 9 


1108 


508 


C-MM4 


06:41:10.0 


+09:29:20 


18x15 


-11 


426 ± 8 


951 


541 


C-MM5 


06:41:10.2 


+09:29:36 


21x11 


135 


261 ± 9 


498 


529 


C-MM6 


06:41:14.8 


+09:27:10 


12x12 




183 ± 6 


218 


100 


C-MM7 


06:41:12.7 


+09:27:01 


15x13 


126 


45 + 6 


73 


68 


C-MM8 


06:41:15.2 


+09:27:31 


14x11 


23 


35 + 6 


45 


89 


C-MM9 


06:41:15.1 


+09:29:36 


21x11 


-1 


94 ± 10 


179 


225 


C-MM10 


06:41:08.5 


+09:29:47 


27x12 


-24 


140 ± 8 


375 


366 


C-MM11 


06:41:08.2 


+09:30:00 


16x11 


103 


110 ± 7 


160 


288 


C-MM12 


06:41:09.8 


+09:29:51 


14x11 


121 


183 ± 7 


233 


387 


D-MM1 


06:41:05.8 


+09:34:06 


17x13 


-24 


257 ± 8 


469 


211 


D-MM2 


06:41:06.3 


+09:35:56 


23x14 


116 


134 ± 9 


357 


204 


D-MM3 


06:41:04.2 


+09:35:01 


18x12 


-11 


105 ± 9 


187 


201 


D-MM4 


06:41:00.3 


+09:35:26 


26x13 


81 


81 + 9 


226 


159 


D-MM5 


06:40:49.4 


+09:34:29 


19x13 


135 


88 + 8 


180 


93 


D-MM6 


06:40:57.9 


+09:36:24 


20x15 


-39 


86 + 9 


213 


119 


D-MM7 


06:41:11.6 


+09:35:32 


28x13 


126 


60 + 8 


180 


116 


D-MM8 


06:41:02.9 


+09:34:31 


15x11 


33 


41 + 9 


56 


136 


D-MM9 


06:41:07.0 


+09:33:31 


31x12 


139 


77 + 7 


237 


137 


D-MM10 


06:41:07.7 


+09:34:18 


16x11 


12 


36 + 8 


52 


136 


D-MM11 


06:40:47.9 


+09:34:42 


19x15 


109 


62 + 8 


146 


90 


D-MM12 


06:40:51.2 


+09:35:24 


12x11 


51 


49 ± 10 


53 


53 


D-MM13 


06:41:08.5 


+09:35:43 


12x11 


-34 


46 + 8 


50 


137 


D-MM14 


06:40:48.7 


+09:35:57 


12x11 


29 


46 ± 10 


50 


26 


D-MM15 


06:41:04.6 


+09:36:19 


16x12 


102 


62 + 9 


98 


146 



a The C-MM and D-MM numbers correspond to our own labelling, and are consistent with the source numbering of Ward- 
Thompson et al. (2000) and Wolf-Chase et al. (2003) for previously known sources. 

b Undeconvolved FWHM sizes derived from fitting an elliptical Gaussian to the background subtracted maps (i.e., after 
filtering out emission seen on scales larger than ScUp = 24") 

c Position angle (from North to East) of the major axis of the fitted Gaussian ellipse 

d Peak flux density in the background subtracted continuum map and la rms at the position of the source 
e Total integrated flux density under the elliptical Gaussian 
' Background level at the position of the source peak 



In order to search for compact millimeter contin- 
uum sources (hereafter, MMSs) within NGC 2264-C and 
NGC 2264-D, we used the systematic method developed 
by Motte, Schilke, & Lis (2003). This method is based on 
a combination of the multi-resolution wavelet algorithm of 
Starck et al. (1995) and the Gaussian fitting procedure of 
Stutzki & Glisten (1990), Gaussclump (kindly provided by 
C. Kramer). The original image is first decomposed into 
two wavelet "views" of the field on small and large spatial 
scales, respectively. The difficult step at this stage is to se- 
lect a proper limiting scale, called ScIJp, beyond which the 
detected dust emission is deemed to arise from the ambi- 
ent cloud rather than from prestellar and/or protostellar 
sources. Here, we choose ScUp by analogy with the previ- 
ous detailed study of the nearby p Ophiuchi protocluster 
by Motte, Andre, & Neri (1998 - hereafter MAN98). In 
p Ophiuchi, the physical diameter of dense cores (such as, 
e.g., Oph-A) was found to be ~ 0.1 pc, corresponding to 
an angular diameter ~ 24" at 800 pc. As the mean an- 



gular separation between nearest neighbouring peaks in 
our dust continuum map of NGC 2264-C is also ~ 24", 
we adopt ScUp = 24" in NGC 2264-C. In order to facili- 
tate the comparison of source masses at the same spatial 
scale in the entire NGC 2264 region, we adopt ScUp = 
24" in NGC 2264-D as well, even though the mean an- 
gular separation between nearest neighbouring peaks is 
slightly larger (~ 40") in clump D. The wavelet decompo- 
sition then allows us to filter out all emission structures 
arising from spatial scales larger than ScUp (see Fig. [2p 
and Fig. Eh). 

In a second step, we use Gaussclump (Stutzki & 
Giisten 1990) to identify all Gaussian sources with FWHM 
sizes larger than the 11" beam width and above a given 
peak flux density threshold in the filtered image (see Motte 
et al. 2003 for details). Following this procedure with the 
detection threshold set to 5<r = 35 mJy/ll"-beam, we 
identify a total of 27 MMSs : 12 in NGC 2264-C (see 
Fig. |3>) and 15 in NGC 2264-D (see Fig. EJd), 14 of 
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(a) NGC 2264-C: 1.2mm dust cont. 



NGC 2264-D: 1.2mm dust cont. 
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Fig. 2. a) Dust continuum map of NGC 2264-C at 
1.2 mm. The white open star shows the position of 
IRS1. The outer white contour corresponds to Sp e 2 fc =70 
mjy/ll"-beam (i.e. 10a), equivalent to a column den- 
sity ~ 7xl0 22 cm~ 2 (assuming k=0.005 cm 2 .g _1 and 
Td=15 K), the black contours go from 120 to 270 by 
50 mJy/ll"-beam and from 350 to 450 by 100 mJy/11"- 
beam, while the inner white contours go from 550 to 1050 
by 100mJy/ll"-beam b) Filtered 1.2mm dust continuum 
image of NGC 2264-C, obtained from the wavelet plane 
corresponding to spatial scales smaller than ScUp=24". 
The white crosses mark the central positions of the MMSs 
identified with Gaussclump. The contours go from 35 to 
50 by 15 mJy/ll"-beam, from 50 to 110 by 20 mJy/11"- 
beam and from 140 to 440 by 50 mJy/ll"-beam (See text 
for details.) 



which are new 1.2 mm dust continuum detections. Tabled 
lists the parameters of the 27 detected Gaussian MMSs. 
Although three more peaks lie just above 35 mJy/beam 
in the filtered map of NGC 2264-D (FigEb), they are 
not identified as real MMSs since their FWHM sizes are 
smaller than the beam width. 
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Fig. 3. Same as FigElfor NGC 2264-D. a) The white open 
star shows the position of IRS2. The outer white contour 
is the same as in Fig. |2t, while the other black contours 
go from 110 to 270 by 40 mJy/ll"-beam and from 350 
to 425 by 75 mJy/ll"-beam. b) The contours go from 35 
to 95 by 10 mJy/ll"-beam and from 120 to 200 by 20 
mJy/ll"-beam (See text for details.) 



3.2. Properties of the NGC 2264 sources 

Assuming optically thin 1.2 mm dust continuum emis- 
sion, the measured flux densities provide direct estimates 
of the masses and column densities of the MMSs (see, 
e.g., MAN98 for details). We adopt a dust mass opacity 
Kdust = 0.005 cm 2 .g _1 at 1.2 mm (cf. MAN98) and a dust 
temperature Td = 15 K (Ward-Thompson et al. 2000) 
for all sources. The distance to NGC 2264 is taken to be 
d = 800 pc and the mean molecular weight fi = 2.33. The 
derived masses, column densities, and volume densities are 
listed in Tabled for both NGC 2264-C and NGC 2264-D. 
Note that the masses of the compact MMSs account for 
only ~ 10% of the total gas mass of the clumps, the latter 
being dominated by the large-scale background. 
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Table 2. Properties of the millimeter sources detected in NGC 2264-C and NGC 2264-D 



Source 


Dec.FWHM a 


M L2 ' 


Density c 






a Un J 


M ■ 9 


a h 




(PC) 


(M Q ) 


(cm" 3 ) 


(10 22 cm- 2 ) 


(10 22 cm- 2 ) 


(kirns" 1 ) 


(M Q ) 




C-MM1 


0.035xUR 


13.1 


1.8xlO b 


26 


19 


0.38 


3.2 


0.2 


C-MM2 


0.076x0.028 


16.0 


6.7xl0 5 


18 


30 


0.51 


8.4 


0.5 


C-MM3 


0.057x0.028 


40.9 


2.7xl0 6 


57 


51 


1.06 


31.3 


0.8 


C-MM4 


0.057x0.041 


35.1 


1.3xl0 6 


43 


54 


0.89 


26.7 


0.8 


C-MM5 


0.072xUR 


18.4 


8.4xl0 5 


26 


53 


0.93 


27.1 


1.5 


C-MM6 


unresolved 


8.1 


>1.5xl0 6 


18 


10 


0.34 


2.3 


0.3 


C-MM7 


0.041x0.028 


2.7 


2.9xl0 5 


5 


7 


0.59 


8.2 


3.0 


C-MM8 


0.035xUR 


1.7 


2.3xl0 5 


4 


9 


0.47 


4.8 


2.8 


C-MM9 


0.072xUR 


6.6 


3.0xl0 5 


9 


23 


0.55 


9.5 


1.4 


C-MM10 


0.099xUR 


13.9 


3.9xl0 5 


14 


37 


0.76 


21.2 


1.5 


C-MM11 


0.046 xUR 


5.9 


5.3xl0 5 


11 


29 


0.68 


11.6 


2.0 


C-MM12 


0.035xUR 


8.6 


1.2xl0 6 


18 


39 


0.76 


12.6 


1.5 


D-MM1 


0.052x0.028 


17.3 


1.3xlO b 


26 


21 


0.59 


9.3 


0.5 


D-MM2 


0.081x0.035 


13.2 


3.6 xlO 5 


13 


20 


0.51 


9.7 


0.7 


D-MM3 


0.057xUR 


6.9 


4.5xl0 5 


11 


20 


0.26 


1.9 


0.3 


D-MM4 


0.094x0.028 


8.4 


2.6 xlO 5 


8 


16 


0.81 


23.5 


2.8 


D-MM5 


0.062x0.028 


6.7 


3.8xl0 5 


9 


9 


0.42 


5.1 


0.8 


D-MM6 


0.067x0.041 


7.9 


2.3xl0 5 


9 


12 


0.55 


11.0 


1.4 


D-MM7 


0.103x0.028 


6.7 


1.8 xlO 5 


6 


12 


0.59 


13.0 


1.9 


D-MM8 


0.041 xUR 


2.1 


2.2xl0 5 


4 


14 


0.34 


2.7 


1.3 


D-MM9 


0.116xUR 


8.9 


2.0xl0 5 


8 


14 


0.55 


12.0 


1.3 


D-MM10 


0.046xUR 


1.9 


1.5xl0 5 


4 


14 


0.72 


13.0 


6.8 


D-MM11 


0.062x0.041 


5.4 


1.7xl0 5 


6 


9 


0.30 


3.2 


0.6 


D-MM12 


unresolved 


1.9 


>3.6xl0 5 


5 


5 


0.38 


2.8 


1.5 


D-MM13 


unresolved 


1.8 


>3.4xl0 5 


5 


14 


0.89 


15.5 


8.6 


D-MM14 


unresolved 


1.8 


>3.4xl0 5 


5 


3 








D-MM15 


0.046 xUR 


3.6 


3.2xl0 5 


6 


15 


0.64 


10.2 


2.8 



a Deconvolved major and minor FWHM continuum sizes; UR means unresolved 

6 Mass estimated from the total 1.2 mm continuum flux density of the corresponding elliptical Gaussian source. For a 2D 
Gaussian source, ~ 94% of this mass is contained within an area twice the size of the FWHM ellipse. Typical uncertainty is a 
factor 2 (on either side) due to the uncertain values of the dust mass opacity and dust temperature 

c Mean H2 volume density derived from the mass of col. [3] assuming a spherical source of effective diameter twice the 
geometrical mean of the deconvolved FWHM sizes in col. [2]. Typical uncertainty is a factor ^ 2 as for Mi. 2. The mean density 
within the FWHM volume is four times larger. 

Peak column density estimated from the background-subtracted map. Typical uncertainty is a factor ^ 2 as for Mi. 2 

e Background column density at the position of the source peak assuming the same dust properties as for the MMSs (see 
text). 

f Velocity dispersion (i.e. rms) along the line of sight derived from a Gaussian hyperfine fit to the N2H + (l-0) multiplet. The 

typical fit uncertainty is Ao"(ine=0.04 km.s" 1 

2 

9 Virial mass calculated as M„, r = 3R assuming a p <x r~ density profile. Here, the radius R is equal to twice the 
geometrical mean of the deconvolved HWHM radii. Typical relative uncertainty is less than 30% 
h Virial parameter defined as a V i r — M„ r /Mi.2. Typical uncertainty is a factor £ 2, dominated by the uncertainty on Mi. 2 



Given the Gaussian fit uncertainties, a source is con- 
sidered as unresolved if its undeconvolved FWHM angular 
size is smaller than 13". Four sources are found to be un- 
resolved (with both the major and minor FWHM sizes 
smaller than 13"), twelve are only partly resolved (with 
only the minor FWHM smaller than 13"), and eleven are 
fully resolved (with both the major and minor FWHMs 
larger than 13"). The MMSs have deconvolved FWHM 
physical sizes ranging from < 0.03 pc to 0.12 pc, masses 
ranging from ~ 2 to ~ 41 M Q , and volume densities rang- 
ing from ~ 1.5xl0 5 to 3xl0 6 cm" 3 (we can only esti- 
mate lower limits to the volume densities of the unre- 



solved MMSs). Our 5(7 detection threshold corresponds 
to a mass of ~ 1.2 Mq. A detailed comparison of the 
global properties of NGC 2264-C and NGC 2264-D is pre- 
sented in § 8 and Table El below. Here, we simply note that 
NGC 2264-C appears to be more centrally concentrated 
than NGC 2264-D in Fig. 1. 

4. Molecular line mapping results 

4.1. N2H + line widths and velocity dispersion 

Maps in low optical depth tracers of dense gas such as 
N2H + can provide constraints on the velocity field of pro- 
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Table 3. Mean N2H + (l-0) line properties of clumps C 
and D 



NGC 2264-C: N 2 H + (1-Q) 



Clump 


V 

V sys 


< &Mne > 




c 


O-iD 


d 




(kirns -1 ) 


(km.s" 1 ) 


(km.s 


- 1 ) 


(km.s 


- 1 ) 


NGC 2264-C 


7.5 ± 0.2 


1.1 ± 0.1 


0.7 ± 


0.1 


1.3 ± 


0.3 


NGC 2264-D 


5.5 ± 0.2 


1.0 ± 0.1 


0.8 ± 


0.2 


1.4 ± 


0.3 



a Mean systemic velocity and rms dispersion around that 
value 

b Mean N2H + (l-0) line-of-sight velocity dispersion averaged 
over each clump and rms dispersion around that value 

c Standard deviation of the distribution of source centroid 
velocities measured in each clump. The error bar has been esti- 
mated as 



assuming that the source sample is drawn 

from a larger population whose velocity distribution follows 
Gaussian statistics. 

3D velocity dispersion of the MMSs calculated from am 
assuming isotropic motions. The error bar has been scaled from 
that estimated for o\d 



toclusters projected onto the line of sight (cf. Belloche 
et al. 2001). The N2H + molecule, which does not de- 
plete up to fairly high densities (e.g. Bergin & Langer 
1997, Belloche & Andre 2004), is a particularly interest- 
ing tracer in this respect. Figure Q] and Fig. [S] show that 
the N 2 H+(l-0) integrated intensity maps of NGC 2264- 
C and NGC 2264-D trace essentially the same structures 
as the 1.2 mm dust continuum maps when the latter are 
smoothed to the same 27" angular resolution. 

We used the HFS (HyperFine Structure) fitting rou- 
tine of the CLASS reduction package from IRAM to fit 
all 7 hyperfine components of the N2H + (l-0) multiplet 
at each mapped position. This provided, for each spec- 
trum, estimates of the centroid velocity, Vlsr, and of the 
FWHM linewidth, AV. 

The N2H+ linewidth measured toward each MMS is 
listed as a velocity dispersion au ne — AV/ ^8 ln(2) in 
Table [3 On average, we find <Ju ne ~ 0.6-0.7 km.s -1 in 
both NGC 2264-C and NGC 2264-D. From these a Une es- 
timates, we can calculate a virial mass, and a virial mass 
ratio (cf. Bertoldi & McKee 1992), a mr = M vir /Mi. 2 , 
for each MMS (cf. Table EJ. Note that the N 2 H+(l-0) 
observations probe material on 27" scale, while the typi- 
cal FWHM size of the dust continuum sources is only ~ 
15". However, since both M V i r and Mi. 2 are derived for 
a diameter ~ twice the FWHM size, our method of es- 
timating a v i r should be reliable. Given the observational 
uncertainties on M„j r and Mi. 2, most MMSs are consis- 
tent with virial equilibrium (i.e., a v i r ~ 1). This is dis- 
cussed further in § 8. The measured values of au ne are 
to be compared with the thermal sound speed, c s , which 
is ~ 0.23 km.s" 1 for molecular gas at 15 K. All MMSs 
have (June > c s , indicating that supersonic turbulence still 
dominates over thermal broadening down to spatial scales 
< 0.1 pc in the NGC 2264 region 1 . The typical Mach 




a(J2000) 

(b) NGC 1.2mm smoothed 




■41 80' 



a(J30aO) 



Fig. 4. a) Map of the N 2 H + (l-0) integrated line inten- 
sity from to 25 km.s" 1 in NGC 2264-C. First contour 
at 2 K.km.s" 1 ; other contours go from 4 to 36 K.km.s" 1 
by 4 K.km.s" 1 . b)Millimeter dust continuum map of 
NGC 2264-C smoothed to the same 27" (HPBW) angular 
resolution as the N 2 H+(l-0) map. Contour levels go from 
0.5 to 1.5 Jy/27"-beam by 0.25 Jy/27"-beam and from 
1.5 to 3 Jy/27"-beam by 0.375 Jy/27"-beam. The white 
crosses mark the positions of the millimeter continuum 
sources identified in 8 3. 



number M = y/3— of the gas within the MMSs is larger 
than 5. On the larger ~ 0.5 — 1 pc spatial scale of the 
NGC 2264 clumps, we estimate the Mach number to be 
M £ 7, based on the linewidth of the mean N2H + (l-0) 
spectrum averaged over the whole extent of NGC 2264-C 
and NGC 2264-D, respectively (cf. Table©. We conclude 
that NGC 2264 is characterized by supersonic motions 
down to the spatial scale of the MMSs. 
We can also estimate the dispersion, ct\d, of the line-of- 



1 In some cases, gravitational collapse may contribute to the 
linewidth and the value of M vir quoted in Table |5] may be 



correspondingly overestimated. The radiative transfer models 
we present in § 6 and § 7 below suggest that infall contributes 
~ 20-50% of the linewidth for C-MM3 and D-MM1. 
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NGC 2264-D: N 2 H + (1 -0) 




"41 1 00 

ff(J200Q) 

NGC 2264-D: L2mm smoothed 




oo" 

«(J2000) 



Fig. 5. Same as FigEJfor NGC 2264-D. a) Contour levels 
go from 2 to 20 K.km.s -1 by 2 K.km.s -1 . b) Contour levels 
go from 0.25 to 1.3 Jy/27"-beam by 0.15 Jy/27"-beam. 



sight velocities of the MMSs relatively to one another in 
each clump. Assuming isotropic motions, we can infer the 
3D velocity dispersion of the MMSs within each clump, 
(J id (see Table |3J). Based on these velocity dis- 
persion estimates, the crossing times of the MMSs across 
their parent clumps are calculated to be t cross ~ 6.0 x 
10 5 yr in Clump C and ~ 6.3 x 10 5 yr in Clump D. 

4.2. A velocity discontinuity at the center of 
NGC 2264-C 

Our mapping in low-optical depth transitions (e.g. 
N 2 H+(l-0) and H 13 CO+(1-0)) revealed a remarkable ve- 
locity discontinuity in the innermost part of NGC 2264- 
C, near the position of the continuum source C-MM3. 
This is illustrated in Fig. H3 which shows N2H + (l-0) 
and H 13 CO + (1-0) position- velocity (PV) diagrams ob- 
served along four different axes. In the PV diagrams taken 
along the East- West (EW) axis going through the con- 
tinuum sources C-MM2, C-MM3, C-MM4 (top two pan- 
els of Fig. Hal, one can clearly see two distinct veloc- 



ity components associated with C-MM2 and C-MM4 at 
~ 6.5 km.s -1 and ~ 8.5 km.s -1 , respectively, both ex- 
tending over more than 50" (i.e., ~ 0.2 pc). These two 
velocity components overlap at the position of C-MM3, 
forming a sharp velocity discontinuity ~ 2 km.s -1 in am- 
plitude. The other panels of Fig. [5] show PV diagrams 
taken along North-South axes at three different right- 
ascension positions, i.e., passing through C-MM2 (labelled 
NS-2 axis), C-MM3 (NS-3 axis), and C-MM4 (NS-4 axis), 
respectively. Along the NS-4 axis, only the higher veloc- 
ity component at ~ 8.5 km.s - can be seen, while along 
the NS-2 axis only the lower velocity component at ~ 6.5 
km.s -1 is visible. Both velocity components are visible in 
the PV diagram along the NS-3 axis. 

We stress that the velocity feature seen in Fig. can- 
not be explained by rotation. Indeed, both velocity com- 
ponents are strong at the center of the system (i.e. ~ C- 
q MM3) with little emission at intermediate velocities, while 
§ the opposite trend would be expected in the case of rota- 
tion. Moreover, the rotational curve expected from differ- 
aential rotation is characterized by a continuous " S" shape 
^.(cf. Belloche et al. 2002) rather than a sharp velocity dis- 
continuity as observed here. 

CO 

i = Finally, we note that no similar velocity discontinuity 
1 o.5g exists at the center of the NGC 2264-D clump. 



4.3. HCO + (3-2) signatures of infall and outflow 

Our extensive HCO+(3-2) maps of both NGC 2264 clumps 
are shown in Fig.0 in the form of spectra overlaid on the 
1.2 mm dust continuum images. The HCO + (3-2) transi- 
tion, which is optically thick and often self-absorbed in 
dense cores, is a good tracer of inward/outward motions 
(e.g. Evans 1999) when associated with an optically thin 
tracer such as H 13 CO+(3-2). A double-peaked HCO+(3-2) 
spectrum with a blue peak stronger than the red peak is 
usually taken to be a diagnostic of infall motions, provid- 
ing the corresponding optically thin H 13 CO + (3-2) spec- 
trum peaks in the dip of the HCO + (3-2) line profile. 
Conversely, an HCO + (3-2) spectrum skewed to the red, 
with, e.g., a red peak stronger than the blue peak, diag- 
noses outflow motions. These diagnostics are valid only 
for centrally-condensed sources (such as dense cores) as 
they assume that the line excitation temperature increases 
toward source center. In NGC 2264-C and NGC 2264- 
D, a variety of HCO + (3-2) line shapes are observed (cf. 
Fig. UJ , ranging from clear 'blue' infall profiles toward C- 
MM3 and D-MM1, to typical 'red' outflow profiles toward 
C-MM4 and D-MM3. Examples of unclear, mixed pro- 
files are observed toward C-MM2 or D-MM2. This variety 
of HCO + (3-2) line shapes suggests that the velocity field 
within the NGC 2264 clumps is complex and possibly re- 
sults from a combination of small-scale inflows/outflows 
directly associated with the MMSs and larger-scale sys- 
tematic motions in the entire clumps. This will be dis- 
cussed in 8 6 and § 7 below. 
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Fig. 7. a) Map of HCO + (3-2) spectra observed toward 
NGC 2264-C, overlaid on the 1.2 mm dust continuum im- 
age (grey scale); the (0,0) position corresponds to C-MM3. 
b) Same as (a) for NGC 2264-D; the (0,0) position corre- 
sponds to D-MM3. 



5. Nature of the embedded sources of NGC 2264 

We have classified the extracted MMSs as "protostellar" 
or "prestellar" based on the presence or absence of an 
outflow, jet or embedded YSO signature (see Tableland 
Tabic 0. For this purpose, we use existing outflow/jet 
studies of the NGC 2264 region (Margulis et al. 1988; 
Schreyer et al. 1997, 2003; Wang et al. 2002), as well as the 
extensive 2MASS and MSX near-/mid-IR surveys. Our 
classification of the NGC 2264-D sources relies entirely on 
comparison with the 2MASS and MSX surveys as no sys- 
tematic jet survey is available for clump D. By default, 
MMSs with no IR counterparts in NGC 2264-D are ten- 



tatively classified as "prestellar" , even though we cannot 
rule out that some may be Class protostellar sources. 

We find that up to ~ 70% of the MMSs in clump C 
and at least 25% of the MMSs in clump D are protostel- 
lar in nature. Although we generally lack good bolometric 
luminosity estimates for the protostellar MMSs and thus 
cannot prove that their submillimeter to bolometric lu- 
minosity ratios match the definition of Class protostars 
(Andre et al. 1993), the fact they exhibit strong millimeter 
continuum emission suggests that they all are good can- 
didate Class objects. We note that the Class nature of 
one of the protostellar MMSs, D-MM1, has been reliably 
established on the basis of SCUBA imaging and HIRES 
processing of the IRAS data (Wolf-Chase et al. 2003). 

The apparent difference in protostellar fraction be- 
tween clump C and clump D may be partly due to the 
lack of outflow /jet survey in clump D. As an illustration, 
C-MM1 and C-MM2 in clump C are both classified as 
protostars because they show evidence of jets, while they 
are not associated with any near-/mid-IR sources. The 
high fraction of protostars found in NGC 2264-C contrasts 
with the significantly lower ~ 40% fraction of protostellar 
sources observed by MAN98 in p-Oph. It is also notewor- 
thy that the protostellar sources of NGC 2264-C are all 
concentrated in the dense, inner ridge visible in Fig. 
while the prestellar sources are all located in the outskirts 
of clump C. 

When comparing our present observations of 
NGC 2264 with the study of MAN98 in p Oph, one 
should keep in mind that there is a difference of a factor 
of ~ 5 in distance, hence in effective spatial resolution. 
Based on the 1.2 mm continuum mosaic of MAN98, if 
p-Oph were at the same distance as NGC 2264 (800 pc), 
our source extraction method would only have detected 
four sources (above the same detection threshold as in 
NGC 2264), corresponding to the prominent dense cores 
Oph A, Oph B2, Oph C, and Oph F (cf. Loren, Wootten, 
& Wilking 1990). We conclude the MMSs identified here 
in NGC 2264 resemble more the DCO + 'dense cores' of 
p Oph than the compact starless condensations found 
by MAN98 (called 'starless clumps' therein). This may 
partly explain why we find a larger fraction of protostellar 
sources in NGC 2264: each dense core in p-Oph harbors 
at least one protostar and several starless condensations; 
observed from a distance of 800 pc, such a core would 
likely be classified as 'protostellar', even though the 
majority of its small-scale condensations actually are 
'prestellar'. 

On the other hand, millimeter continuum interferom- 
eter observations by Nakano et al. (2003) and Schreyer et 
al. (2003) show that the objects called C-MM4, C-MM5, 
and C-MM12 here remain unfragmented, single sources at 
a spatial resolution of ~ 3000 AU. This may suggest that 
the MMSs detected with the 30 m telescope in NGC 2264 
are more centrally concentrated than the p-Oph dense 
cores (see § 8 for further comparisons). 
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Table 4. Nature of the millimeter continuum sources iden- 
tified in NGC 2264-C 



Source MSX/2MASS° 


H 2 jet" 


CS fiow c Nature" 


C-MM1 


Y 


pro 


C-MM2 


Y 


pro 


C-MM3 3 


Y 


pro 


C-MM4 1 


? 


N pro 


C-MM5 2 


? 


N pro 


C-MM6 


N 


pre 


C-MM7 


N 


pre 


C-MM8 


N 


pre 


C-MM9 


N 


pre 


C-MM10 1 


N 


Y pro 


C-MM11 2 


N 


pro 


C-MM12 2 


? 


Y pro 


a Number of 2MASS or MSX infrared 


sources lying within 


an 11" beam centered on the MMS position 


6 Detection of a shocked H2 


jet by Wang et al. (2002) 


c Detection of a CS outflow by Schreyer et al. (2003) 


d 'pro' stands for protostellar, 'pre' for prestellar 


Table 5. Nature of the millimeter continuum sources iden- 


tified in NGC 2264-D 






Source MSX/2MASS a 


Nature" 


D-MM1 


1 


pro 


D-MM2 





pre 


D-MM3 


1 


pro 


D-MM4 





pre 


D-MM5 





pre 


D-MM6 





pre 


D-MM7 


1 


pro 


D-MM8 





pre 


D-MM9 





pre 


D-MM10 





pre 


D-MM11 





pre 


D-MM12 





pre 


D-MM13 





pre 


D-MM14 


1 


pro 


D-MM15 





pre 



a Number of 2M ASS or MSX infrared sources lying within 
an 11" beam centered on the MMS position 
b 'pro' stands for protostellar, 'pre' for prestellar 



The bright infrared source IRS1 lies close to, but is 
clearly offset from, C-MM5 (Nakano et al. 2002; Schreyer 
et al. 2003). Based on its slightly rising spectral en- 
ergy distribution between 12 /im and 100 /im, Margulis 
et al. (1989) classified IRS1 as a Class I object with 
Lfc o z=2300 Lq. The detections of VLA radio continuum 
emission (Schwartz et al. 1985), as well as H2O and 
methanol masers, at the position of IRS1 are suggestive 
of a - 9.5 Mq YSO on the ZAMS (cf. Thompson et al. 
1998). Although there is apparently no disk nor outflow 
directly associated with it (Schreyer et al. 2003), IRS1 is 
thus clearly a relatively massive, embbeded young star. 
The other infrared source detected by IRAS in the region 



discussed in this paper is IRS 2, which lies in clump D, 
just outside the densest region (see Fig. 1 and Fig. 13). 
It is significantly less luminous than IRS1 (Lf, o j = 150 L Q 
- Margulis et al. 1989) and is not closely associated with 
any millimeter continuum source. It has also been classi- 
fied as a relatively young (Class I) object by Margulis et 
al. (1989). 

6. Radiative transfer modelling of NGC 2264-C 

Here, we compare our molecular line observations of 
NGC 2264-C with radiative transfer calculations per- 
formed with the same Monte-Carlo radiative transfer code 
as used by Belloche et al. (2002). This code is divided into 
two parts. The first part calculates the non-LTE level pop- 
ulations based on a ID Monte-Carlo method (Bernes 1978, 
1979). The second part, called MAPYSO (Blinder 1997), 
integrates the radiative transfer equation along each line 
of sight and convolves the resulting data cube with the 
beam of the telescope. MAPYSO can be used in either ID 
or 2D geometry. The 2D mode was used for the present 
modelling of NGC 2264-C. The inputs to the code are the 
density and abundance profiles, the kinetic temperature 
profile, the 2D velocity field, and the non-thermal velocity 
dispersion which, for simplicity, is assumed to be uniform 
over the simulated region. 

The averaged radial density profile of NGC 2264-C can 
be inferred from the corresponding 1.2 mm dust contin- 
uum intensity profile I(r), where r is the spherical radius. 
The best power-law fit to the observed intensity profile is 
found to be I(r) oc r~ 4 . Assuming isothermal dust emis- 
sion in the Rayleigh- Jeans regime, this translates into an 
average spherical density profile p(i) oc r -1 ' 4 (cf. Motte & 
Andre 2001), which has been normalized so as to yield a 
total mass of 1650 M Q within a radius of 0.4 pc, as ob- 
served (cf. Table El below) . Since the projected center of 
gravity of the NGC 2264-C clump as estimated from our 
dust continuum map is located at the position of C-MM3, 
we have centered our models on C-MM3 rather than C- 
MM4 as adopted by Williams & Garland (2002). 

The gas kinetic temperature profile is taken to be 
T k (r) = ((T (^)-°- 4 ) 4 + T 4 ) - 25 with T o =90 K, r =100 
AU and Ti = 15 K. The value of To has been adjusted so 
as to match the intensity of the observed spectra. This 
temperature profile is consistent with heating from a lu- 
minous central source (Lboi ~ 80 L©), radiating in an 
optically thin medium and embedded in a cold 15 K back- 
ground at large radii. 

The CS and HCO + molecules are known to be de- 
pleted onto dust grains above a density ~ 10 5 cm~ 3 
in starless cores (e.g. Tafalla et al. 2002). The spher- 
ical abundance profiles are thus fixed to standard val- 
ues, i.e. [CS]/[H 2 ] = lxl0- 9 (i.e. Tafalla et al. 2002) and 
[HCO+]/[H 2 ] = 1.5xlO- 9 (i.e. Bergin et al. 1997), in the 
low-density outer regions and are assumed to scale as 
n H° 2 4 w ith density (cf. Bacmann et al. 2002). We assume 
the decrease in abundance with density to stop where 
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(b) 



Fig. 8. a) Schematic view of our proposed model of 
NGC 2264-C as seen by an observer in the plane of the 
sky. b) Side view of the plane containing the z axis and 
the line of sight: a finite-sized, cylindrical clump collapses 
toward its center, coinciding with C-MM3. The long axis 
(i.e. z axis) of the cylinder makes an angle of 45 degrees 
with the line of sight. 

the kinetic temperature rises above 20 K in the inner re- 
gion (see Fig- El - Still further in, where the gas becomes 
warmer than ~ 50 K, one may expects the abundances to 
increase as a result of grain mantle evaporation. However, 
Tfc reaches 50 K only at very small radii (< 1000 AU cor- 
responding to < 1") in our model, and this small inner 
region has no influence on spectra observed at > 10" an- 
gular resolution. We also adopt constant, standard values 
of 22 and 75 for the [CS]/[C 34 S] and [HCO+]/[H 13 CO+] 
isotopic ratios, respectively. 

The non-thermal velocity dispersion is set to cfnt— 0.47 
km.s -1 , i.e., slightly less than the N2H + (l-0) linewidths 
measured in NGC 2264-C (cf. § 14. lj) . which accounts 
for the contribution of collapse motions to the linewidth. 
Finally, the infall velocity profile is adjusted so as to fit 
the observed spectra. 

In order to account for the velocity discontinuity ob- 
served in the central part of NGC 2264-C, we model the 
clump as a cylindrical filament centered on C-MM3 (see 
Fig.lHJ). The system is assumed to be collapsing toward its 
center of mass (i.e. C-MM3), so that the two sides of the 
filament (C-MM2 and C-MM4) are moving toward each 
other. The long axis of the filament is taken to be inclined 
by 45° to the line-of-sight (see Fig. |Ht>) ■ When probing the 
kinematics of such structure, one would indeed observe a 



Radius (arcsec) Radius (arcsec] 

1 10 1 10 2 1 10 1 10 1 
^10 a ! ■ i 
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Radius (AU) Z axis (AU) 



Fig. 9. Profiles of the input parameters adopted in our 2D 
radiative transfer model of NGC 2264-C : Density scaling 
as n oc r~ 14 ; depletion of the CS abundance by a factor 
of - 5 in the inner region (the C 34 S, H 13 CO+, and HCO+ 
abundances are assumed to scale accordingly); gas kinetic 
temperature scaling as T^ cx r~ 04 at small radii and fixed 
to Tfe = 15 K in the outer region; velocity profile along 
the z-axis of the cylinder (solid lines), with the locations 
of iVi n f and ±Ri n f marked by dash-dotted lines. 

sharp velocity discontinuity near the center of the fila- 
ment. The velocity profile is taken to be only a function 
of z (position along the long axis of the cylinder), with 
v(z)=Vj n /=1.3 km.s -1 for z<0 and v(z)=-Vj n / for z>0 
in the inner \z\ < 35000 AU (or 45") region, and sharp, 
symmetrical decreases of \v(z)\ with \z \ outside the central 
region (cf. Fig.EJ). 

Figure compares the synthetic spectra calculated 
under these assumptions with the spectra observed along 
the EW axis in the CS(3-2), CS(5-4), C 34 S(2-1), C 34 S(3- 
2), C 34 S(5-4), HCO+(3-2), H 13 CO+(3-2), and H 13 CO+(l- 
0) transitions (where the (0",0") position corresponds to 
C-MM3). The agreement between the simulations and the 
observations is very encouraging given the simplicity of the 
model. The linewidths, peaks, and dips of the observed 
spectra are generally well reproduced. Furthermore, our 
simple cylindrical model manages to reproduce the ob- 
served reversal in line asymmetry, from blue-skewed spec- 
tra near C-MM3 and on the Eastern side of C-MM3 to 
red-skewed spectra on the Western side of C-MM3 (cf. 
Fig. I10|) , without including any additional outflow veloc- 
ity component near C-MM4. On the Western side of C- 
MM3, however, the intensity of the synthetic spectra is not 
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Isocontours of 




Fig. 11. Schematic view of the collapsing NGC 2264-C 
filament explaining the reversal of self-absorbed asymme- 
try observed in optically thick transitions (see § 6). The 
positions of C-MM2, C-MM3, C-MM4 are marked. The 
white dashed circles represent isocontours of increasing ex- 
citation temperature toward C-MM3. Three lines of sight 
(l.o.s.) are shown and the point where the optical depth 
r = 1 along each of them indicated. 



strong enough, especially in the case of the optically thick 
HCO + (3-2) and CS(5-4) transitions. This is because our 
ID model underestimates both the actual density and the 
actual kinetic temperature, hence the excitation tempera- 
ture, near the luminous IR source IRS1 and the millimeter 
sources C-MM4 and C-MM5 (which are all offset from the 
center of gravity of the system at C-MM3). 
We note that both the observed and synthetic H 13 CO + (l- 
0) spectra are asymmetric and even double-peaked near 
C-MM3. This is somewhat reminiscent of the asymmetry 
observed in HCO + (3-2), although there is a fundamen- 
tal difference: While the HCO + (3-2) spectra are optically 
thick and self-absorbed (with an optical depth r ~ 45 in 
the central velocity channel of the model), the H 13 CO + (l- 
0) transition is essentially optically thin (with r ~ 0.6 in 
the central velocity channel of the model), so that the 
asymmetric, double-peaked H 13 CO + (1-0) profiles reflect 
the presence of two velocity components along the line- 
of-sight (see also the position-velocity diagrams shown in 
Fig. El rather than self-absorption effects. 

Figure ITT1 provides a schematic explanation of the line 
reversal for optically thick spectra. Three lines of sight 
(l.o.s. 1, l.o.s. 2, and l.o.s. 3) are represented which inter- 
sect different isocontours of excitation temperature. Due 
to the elongated, cylindrical structure of the model, the 
blueshifted material probed along l.o.s. 3 has much lower 
excitation and/or column density than the blueshifted 
material along l.o.s. 2 or l.o.s. 1. This explains why the in- 
tensity of the blue peak of the self-absorbed lines is strong 
on the Eastern side of C-MM3 and decreases rapidly on 
the Western side (see Fig. By contrast, the two lines 



of sight l.o.s. 2 and l.o.s. 3 include similar amounts of 
excited, redshifted material, which qualitatively explains 
why the red peak remains approximately constant (in 
the model) on the Western side of C-MM3. Of course, 
protostcllar outflows are likely to play an additional 
role in shaping the line profiles observed in this region, 
especially around C-MM4 and IRS 1 (see the broad 
wings of the central CS(3-2) spectrum observed toward 
C-MM4 in Fig. [TOjl . Nevertheless, the above qualitative 
explanation of the asymmetry reversal and the reasonably 
good model fit shown in Fig. HOI suggest that the velocity 
field near C-MM3 is dominated by the global collapse 
of the entire NGC 2264-C elongated clump rather than 
by local inflowing/outflowing motions around individual 
sources. 

In order to estimate rough error bars on our 'best-fit' 
model shown in Fig. 1101 we have performed a set of 
radiative transfer calculations by varying some of the 
model parameters. We find acceptable fits for an inflow 
velocity Vi„/=1.3 km.s -1 ± 0.2 km.s -1 and an inflow 
radius R m/ =35000 AU ± 10000 AU (see Fig. Et for 
definitions of these two parameters). As an illustration, 
Fig. ^] shows synthetic H 13 CO + (1-0) position-velocity 
diagrams for three different models: (a) the "best-fit" 
model, (b) a model similar to the best-fit model but 
with R m/ =20000 AU instead of R m/ =35000 AU , and 
(c) a model for which the velocity field corresponds 
to (differential) rotation about C-MM3 as opposed to 
inflow toward C-MM3. It can be seen that the agreement 
between the best-fit PV diagram (Fig. 112b ) and the 
observed diagram (top right of Fig. EJ) is quite good, 
with two separate velocity components overlapping at 
the central position. By contrast, the model with a 
smaller inflow radius (Fig. 112b ) does not match the 
main features of the observed PV diagram, even if an 
unresolved velocity discontinuity is apparent in the 
center. Likewise, the rotation model does not fit the 
observed PV diagram, producing a continuous velocity 
gradient as opposed to a velocity discontinuity. Based on 
the rotational models we have calculated, we can rule 
out rotation as the origin of the kinematic properties of 
NGC 2264-C. We also tried but failed to reproduce our 
observations with the kinematical model proposed by 
Williams & Garland (2002), namely large-scale spherical 
infall onto an expanding central core coinciding with 
C-MM4. Our higher-resolution observations show that 
the large-scale collapse motions in NGC 2264-C converge 
toward a position much closer to C-MM3 than to C-MM4 
and that these motions are not disrupted on small scales 
by the effect of protostellar outflows around C-MM4 as 
proposed by Williams & Garland (2002). 
Based on our best-fit model (cf Fig. 110(1 , the dynamical 
timescale of the inner part of NGC 2264-C is estimated to 
be tdyn ~ 1.7x 10 5 yr, which corresponds to the character- 
istic time needed by C-MM2 and C-MM4 to collapse to the 
system central position (C-MM3). We can also calculate 
the mass inflow rate toward the center of the cylindrical 
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Fig. 12. Synthetic H 13 CO+(1-0) position-velocity dia- 
grams calculated for three different radiative transfer 
models, a) Diagram corresponding to the 'best-fit' model 
(see text for details), b) Diagram obtained by changing 
Kinf to 20000 AU. c) Diagram obtained for a bulk (differ- 
ential) rotation of the model cloud about an axis centered 
on C-MM3 and perpendicular to the line joining C-MM2, 
C-MM3, and C-MM4. These diagrams have to be com- 
pared with the observed H 13 CO + (1-0) diagram shown at 
the top right of Fig. 



filament M inf = 2 x it R 2 fil x n mean X fj, X m X V in f, 
where Rfu is the radius of the cylinder cross section, 
n me an is the mean number density in the cylinder, /i 
is the mean molecular weight, m is the mass of atomic 
hydrogen, and Vi„/ is the inflow velocity used in our 
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Fig. 13. Radial profiles of the input parameters adopted 
in our 2D radiative transfer model of D-MM1: Density 
scaling as n oc r -1,5 in the inner 10000 AU and n oc r~ 2 
beyond 10000 AU; depletion of the CS abundance by a 
factor of ~ 6 in the inner region (the profiles of the C 34 S, 
HCO + , and H 13 CO + are assumed to have the same form); 
kinetic temperature scaling as T* oc r~ - 4 ; infall veloc- 
ity profile (marked Vi n f) and rotational velocity profile 
(marked V rot and depending only on the cylindrical ra- 
dius from the rotation axis) 



V»n/ = 1.3 km.s , we find 
-l 



model fit. With Rju = 0.2 pc, n„ 



lxlO 5 



(cf Table E5J}, fi = 2.33 

Mi n f ~ 3xl0 -3 Mo.yr -1 . This is an order of mag- 
nitude larger than the mass inflow rate found by 
Williams & Garland (2002), which is not surprising 
since, in our picture, their spherical model was offset 
by ~ 0.15 pc from the true center of mass of NGC 2264-C. 



7. Radiative transfer modelling of D-MM1 

As the mean separation between nearest MMSs is larger 
in NGC 2264-D than in NGC 2264-C, our molecular line 
mapping observations of clump D should be more sensitive 
to the dynamics of individual MMSs. In this section, we 
present an attempt at modelling the strongest millimeter 
source of NGC 2264-D, D-MM1, which was classified as a 
Class protostar by Wolf-Chase et al. (2003) and exhibits 
blue infall profiles in both HCO+(3-2) and CS(3-2) (see, 
e.g., Fig. [J. 

We have modelled D-MM1 as a spherical cloud core 
with both infall and rotational motions using the same 
radiative transfer code (MAPYSO) as for NGC 2264-C 
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(§ 6). Several input parameters of the code are well con- 
strained. We assume the gas kinetic temperature to be well 
coupled to the dust temperature and adopt a dust tem- 
perature profile T d oc r~ 4 with T d ~ 80 K at r=100 AU 
(cf. Motte & Andre 2001 and references therein), consis- 
tent with the observed bolometric luminosity Lb i ~ 100 
Lq (Wolf-Chase et al. 2003). From our 1.2 mm contin- 
uum map, we derive a circularly-averaged radial flux den- 
sity profile which is consistent with a radial density profile 
p oc r -15 in the inner 1 x 10 4 AU radius region, and steeper 
outside with p oc r~ 2 . The model density profile has been 
normalized such that the mass of the cloud core is 10 M Q 
within a radius of 4500 AU as observed. The non-thermal 
velocity dispersion is set to <t^t — 0.51 km.s -1 which 
is slightly less than the observed linewidth (cf Table [2J . 
This is due to the collapse contribution to the linewidth. 
The N2H + (l-0) spectra we observe around D-MM1 in- 
dicate the presence of a relatively strong velocity gra- 
dient ~ 3 km.s^.pc -1 from south-west to north-east 
(P.A.grad ~ +45°). In agreement with our N2H" 1 " data, the 
projection of the rotation axis on the plane of the sky has 
thus been fixed at P.A. rot ~ —45°, and an angular veloc- 
ity f2=3 km.s _1 .pc _1 has been adopted at r = 2 x 10 4 AU 
from the center of D-MM1. 

Taking the above constraints into account, the mag- 
nitudes of the infall and rotational velocities have been 
adjusted so as to match the observations. In the model 
shown in Fig. ED an d Fig. El the infall velocity is such 
as Vi n f oc r~ a5 and fixed to a value of 0.2 km.s" 1 at 
r = 4000 AU. In order to reproduce the spatial sequence of 
self-absorbed HCO + and CS spectra (clearly skewed to the 
blue at negative offsets but almost symmetric at positive 
offsets - see Fig. 1140 , we find that strong differential rota- 
tion is required, with, e.g., solid-body rotation, Y ro t oc r, 
in the inner 2 x 10 4 AU radius region and a sharp decline, 
Y ro t oc r in the outer region. Such a rotational velocity 
profile is reminiscent of the velocity field found by Belloche 
et al. (2002) for the low-mass Class object IRAM 04191 
in Taurus. With these parameters and assumptions (sum- 
marized in Fig. 113(1 . a reasonably good fit is found (see 
Fig. I14f> . although the model fails to reproduce the broad 
wings observed in HCO + (3-2) and CS(3-2), probably due 
to outflowing material. 

Based on this model fit, we can estimate 
the mass infall rate onto D-MM1 as follows: 
Mdmmi — 47rR 2 x uh 2 x p x m x Vj n /. At a radius 
of 10000 AU, V in f ^ 0.1 km.s" 1 and n Ha - 6xl0 5 cm" 3 , 
which yields Mdmmi = l.lxl0~ 4 M .yr _1 . This is 
almost two orders of magnitude larger than the typical 
infall rate measured toward low-mass protostars in 
Taurus (Ohashi 1999, Belloche et al. 2002), and compa- 
rable to the large infall rate estimated for some Class 
protostars in protoclusters such as NGC 1333-IRAS4A 
(cf. Di Francesco et al. 2001, Andre et al. 2004). The 
implications of such a high infall rate are discussed in the 
next section. 
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Fig. 14. CS(3-2), CS(5-4), C 34 S(2-1), C 34 S(3-2) (top), 
and HCO+(3-2), H 13 CO+(3-2), H 13 CO+(1-0) (bottom) 
spectra observed along the direction of maximum velocity 
gradient (P.A. grad ~ +45°) through D-MM1. Synthetic 
spectra corresponding to the "best-fit" 2D radiative trans- 
fer model described in the text are superimposed in light 
grey. Positive offsets are located to the North-East of D- 
MM1, while negative offsets are to the South-West. 

8. Discussion and Conclusions 

8.1. Comparison of three protoclusters: NGC 2264-C, 
NGC 2264-D, and p-Ophiuchi 

The overall dynamical state and stability of a clump may 
influence the formation process of individual stars in its 
interior. We have shown in § 6 that the NGC 2264-C 
clump is in a state of global collapse and evolves on the 
dynamical timescale td yn ~ 1.7xl0 5 yr. This implies 
that NGC 2264-C is not in approximate hydrostatic 
equilibrium. Independently, we can compare the 3D 
velocity dispersion of the MMSs, 03.D, to the 3D velocity 
dispersion expected in virial equilibrium, a v i r _3D- These 
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parameters are given in Table |H| for both NGC 2264-C 
and NGC 2264-D, as well as p-Ophiuchi, the nearest 
(d ~ 150 pc) example of a cluster-forming clump. Table 
also lists the free-fall time of each clump estimated 
from the observed mean volume density. In all three 
cases, the observed velocity dispersion of the MMSs 
appears to be insufficient to balance gravity. Likewise, the 
estimated virial mass M^ is significantly smaller than 
the total gas mass derived from the 1.2 mm continuum 
map for both clumps (cf. Table EJ) , corresponding to a 
virial parameter a V i r — M V i r / Mx.% ~ 0.2, which suggests 
that NGC 2264-C and NGC 2264-D are both out of 
virial equilibrium and highly unstable to collapse. If 
the clumps are gravitationally unstable, we expect the 
free-fall time to provide a reasonably good estimate of the 
characteristic evolutionary timescale. Our independent 
estimate of td yn for NGC 2264-C (see above) supports 
this view since t c i yn is only slightly larger than tff. 
We conclude that the central part of NGC 2264-C is 
close to being in free-fall collapse. NGC 2264-D and 
p-Ophiuchi may also be collapsing. Using HCO + (3-2), 
Williams & Garland (2002) found some evidence for 
large-scale collapse in NGC 2264-D at a characteristic 
velocity v co «. = 0.3 km.s , although more observations 
of optically-thick tracers on large scales would be required 
to confirm their result. Likewise, tentative signatures of 
large-scale infall have been observed in 12 CO and 13 CO 
toward the p Oph clump (e.g. Encrenaz et al. 1975, 
Wilking & Lada 1983). 

We now compare the MMSs of NGC 2264 with the 
DCO + cores of p Oph (Loren et al. 1990) and typical iso- 
lated prestellar cores (Ward-Thompson et al. 1999, Caselli 
et al. 2002). Figure El shows a plot of mean pressure, 
Pcore/k, versus non-thermal velocity dispersion, ct/vt, for 
these various types of cores as well as a typical singular 
isothermal sphere (SIS) model (Shu et al. 1987) and a typ- 
ical turbulent core model (McKee & Tan 2003 - hereafter 
MT03). The mean pressure in each core has been calcu- 
lated from its mean column density derived from millime- 
ter continuum observations, using Eqs. (A5-A6) of MT03 
and assuming a V i r = 1 (which is typical of the NGC 2264 
MMSs - see Table |2J. The non-thermal component of the 
ID velocity dispersion in each core has been estimated 
from the observed N2H" 1 " (101-012) linewidth after sub- 
tracting the thermal broadening expected for Tfc=10 K. 
It can be seen in Fig. ll5l that P CO re/ k increases by three or- 
ders of magnitude and (Jnt increases by one order of mag- 
nitude from the thermally-dominated, isolated prestellar 
cores (cf. Myers 1998) to the most extreme MMSs of 
NGC 2264 (this paper) dominated by supersonic non- 
thermal motions. It can also be noted that the hydro- 
static SIS model provides a reasonably good represen- 
tation of the locations of isolated prestellar cores in the 
Pcore/k vs. ctatt diagram, while the turbulent core model 
of MT03 (see § 8.2 below) provides a much better de- 
scription of the NGC 2264 sources. Interestingly, the low- 
pressure, thermally-dominated cores on the left-hand side 



10 
10 9 

to" - " 1 

E 10 s 

o 

S 1 



□ NGC2264-C 

O NGC2264-D 

A p — Oph Cores 

O Isolated Cores 

• 15K SIS Core 

■ MT03 Core 



L1544 



10 r 



O 



o 

OphA-" 

A 

A 



CMM3 - 
□ 



□ □ 

DM M 1— >0 □ 

□ ■ 

□ 

□ O 

a 



o 



•A °n o 

\ 

OphB2 



(10K) 



N,H4- ( kr 



Fig. 15. Plot of mean pressure versus non-thermal veloc- 
ity dispersion for the millimeter continuum sources iden- 
tified in NGC 2264-C (square symbols) and NGC 2264-D 
(pentagonal symbols), the p Oph DCO + cores Oph A, 
OphB 2, Oph C, Oph F (open triangles - from MAN98 
and Belloche et al. 2001), and the isolated prestellar cores 
L1512, L1544, L183, L1696A, L63, L1155C (open circles - 
from Ward-Thompson et al. 1999 and Caselli et al. 2002). 
The location of a singular isothermal sphere (SIS) model 
with an effective temperature of 15 K and an actual gas 
temperature of 10 K is shown as a filled circle for com- 
parison. The filled square corresponds to the turbulent 
core model of McKee & Tan (2003 - MT03) for a clump 
surface density E c ; = 1 g.cm~ 2 and a final stellar mass 
rn*/ = 10 M Q (see also § 8.2). The vertical dashed line sep- 
arates cores dominated by thermal motions (on the left) 
from cores dominated by supersonic non-thermal motions 
(on the right), assuming a typical gas temperature of 10 K. 



of Fig. 1151 are typically forming one (or two) stars, while 
the high-pressure, "turbulent" cores on the right-hand side 
are associated with protoclusters. A similar distinction be- 
tween the "isolated" and "clustered" modes of star forma- 
tion was noted earlier by Myers (1998) and Jijina et al. 
(1999). 

In the high-pressure, high-cr/vr part of Fig. 1151 it can be 
seen that the p-Oph cores Oph A and Oph B2 are lo- 
cated in the regime of the MMSs of NGC 2264-D. The 
sources of NGC 2264-C have higher mean pressure and 
larger non-thermal velocity dispersions by a factor ~ 10 
and ~ 2, respectively, on average than OphA and OphB2. 
The most extreme_source of NGC 2264-C (C-MM3) has 
a mean pressure P core /k ~ 10 10 Kcm~ 3 reminiscent of 
the pressure found in high-mass star-forming clumps (cf. 
MT03). 

We conclude that the MMSs of NGC 2264, especially 
those of clump C, represent somewhat more pressurized 
and more turbulent conditions for star formation than do 
the p Oph cores. The more extreme conditions found in 
NGC 2264-C are actually reminiscent of those suggested 
by MT03 for the precursors of massive stars. 
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Table 6. Global properties of NGC 2264-C and NGC 2264-D compared to those of p-Oph 



CI. Name 


Cl.Diam. a 


Cl.Mass" 


M - c 


CI. Col. Dens. 11 


Cl.Vol.Dens. c 


t// 
(10 5 yr) 


aso 9 


h 

<?vir_3D 




(pc) 


(M ) 


(Mq) 


(10 22 cm- 2 ) 


(10 4 cm- 3 ) 


(km.s -1 ) 


(km.s -1 ) 


NGC2264-C 


0.8 


1650 


340 


16 


10.7 


1.0 


1.3 


3.7 


NGC2264-D 


0.9 


1310 


310 


11 


6.0 


1.4 


1.4 


3.0 


p-Oph 1 


1.1 


550 




3 


1.4 


2.9 


0.6 


2.1 



a Mean clump diameter measured at the 70 mjy/ll"-beam contour in the 1.2 mm continuum map (cf. Fig. [3^ and Fig- EJ 3 ) 
6 Total gas mass of the clump estimated from the 1.2 mm continuum image prior to spatial filtering 

c Virial mass calculated in the same way as in Table H from the mean N2H + (l-0) line-of-sight velocity dispersion observed in 
the clump (cf. col. [3] of Table |3J 
d Mean column density averaged over each clump 
e Mean density of each clump assuming a spherical volume 
f Mean free-fall time of the clump 

9 Observed 3D velocity dispersion of the MMSs in each clump 
Expected 3D velocity dispersion of the MMSs if they were virialized in the gravitational potential of the clump: ow_3.d = 

' Based on the C 18 O(l-0) results of Wilking & Lada (1983) for the diameter and the mass, and on Belloche et al. (2001) for 
the velocity dispersion 



8.2. Comparison with two scenarios of clustered star 
formation 

The young near-infrared star cluster associated with 
NGC 2264 (Lada et al. 1993) is known to be highly struc- 
tured and exhibits a surface density distribution with mul- 
tiple peaks (see Fig. 4 of Lada & Lada 2003). Two of these 
peaks lie close to, but are clearly offset from, the cluster- 
forming clumps NGC 2264-C and NGC 2264-D studied in 
the present paper (see Fig. lTo|) . This is suggestive of at 
least two different episodes of (clustered) star formation 
in the NGC 2264 region. The detection of a few 2MASS 
and MSX infrared sources toward both NGC 2264-C and 
NGC 2264-D shows that both clumps have already formed 
YSOs. Furthermore, the presence of the massive (B2-type) 
Class I infrared YSO IRS1 within the NGC 2264-C clump 
sets interesting constraints on the process of intermediate- 
to high-mass star formation. 

According to the scenario proposed by Bonnell et 
al. (1998, 2002), the formation of massive stars occurs 
through competitive accretion and stellar mergers in the 
dense inner core of contracting protoclusters when their 
central stellar density exceeds ~ 10 s stars/pc 3 . The grav- 
itational potential well of such cluster cores is dominated 
by the stellar rather than the gas component. We find that 
it is difficult to account for the formation of IRS1 in the 
context of this scenario. Indeed, IRS1 appears to be lo- 
cated in the outer parts of the NGC 2264-C gas clump 
(e.g. Thompson et al. 1998, Nakano et al. 2002), rather 
near the peak of the near-IR source density distribution 
(cf. Fig. I16fl . Figure 1 from Lada et al. (1993) indicates 
a stellar surface density ~ 80 pc~ 2 , which allows us to 
estimate the stellar contribution to the gravitational po- 
tential within the volume of NGC 2264-C. Assuming the 
Kroupa (2001) IMF, for which the mean stellar mass is 
~ 0.4 Mq, the above number surface density translates 
into a stellar mass surface density ~ 30 M Q .pc -2 , i.e., less 



than 2% of the total gas mass surface density. Even allow- 
ing for a possible fraction of highly extinguished stars in 
the protocluster that would not be included in the above 
estimate, we conclude that the gravitational potential of 
NGC 2264-C is dominated by the gas rather than by the 
stars. Thus, IRS1, the most massive infrared YSO in this 
region is unlikely to have formed by collision and coa- 
lescence of lower-mass stars as described by Bonnell et 
al. On the other hand, the global, near free-fall collapse 
of NGC 2264-C may lead to strong dynamical interac- 
tions between MMSs and even protostellar mergers near 
the center of the gas-dominated clump (see below). 

In the alternative scenario proposed by McKee & Tan 
(2002, 2003 - e.g. MT03), high-mass stars form from mas- 
sive, turbulent cores embedded in a high-pressure, envi- 
ronment/clump. Both the cores and the ambient medium 
are assumed to be centrally-concentrated and initially 
virialized, i.e., close to hydrostatic equilibrium. Although 
clumps C and D are indeed massive and turbulent, our ob- 
servations reveal that they are not in hydrostatic equilib- 
rium but rather in a state of global dynamical collapse, at 
least in the case of NGC 2264-C. The primary characteris- 
tic of the MT03 model, however, is that the prestellar pre- 
cursors of massive protostars are high-density, turbulent 
cores, which leads to large, accelerating accretion rates 
during the protostellar phase. Large ( ^ 10~ 4 M Q .yr -1 ) 
accretion rates are indeed necessary to overcome radiation 
pressure and allow high-mass (M* ^ 8M ) stars to form 
by accretion despite the high luminosity of the central 
stars (e.g. Wolfire & Cassinelli 1987). Large-scale collapse 
can help to produce dynamically pressurized, high-density 
cores in the inner part of massive cluster- forming clumps. 
Even though one of the premises of the MT03 model is 
not verified, a detailed comparison of the model with the 
observed properties of NGC 2264-C and its central source, 
C-MM3, is therefore still warranted. 
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Fig. 16. Contour map of the surface density of near- 
IR sources in the young star cluster associated with 
NGC 2264 (from Lada & Lada 2003), overlaid on our 
1.2 mm continuum mosaic of the region (greyscale and 
white contours). The white filled triangles mark the posi- 
tions of IRS1 and IRS2. Note that the maxima of the near- 
IR source surface density distribution are located outside 
the dense gas/dust clumps NGC 2264-C and NGC 2264- 
D, while the projected position of the luminous near-IR 
source IRS1 lies within the 1.2 mm contours of clump C 

The mean column density of the NGC 2264-C clump 
is Njj2 ~ 1.6 x 10 23 cm~ 2 , corresponding to a mean sur- 
face density E c ; ~ 0.6 g.cm~ 2 , typical of that observed in 
most regions of massive star formation and consistent with 
the assumptions of the turbulent core model of MT03. 
With a mass of 40 M© and a line velocity dispersion 
o~Mne — 1.1 km.s -1 , C-MM3 is located near the center 
of gravity of the collapsing NGC 2264-C clump. An up- 
per limit to the mass accretion rate onto the protostellar 
object C-MM3 is provided by the total mass inflow rate 
Minf ~ 3xl0 -3 M Q .yr _1 derived in § 6 for the central 
part of the NGC 2264-C. A lower limit can be estimated 
from the mass accretion rate onto the less extreme object 
D-MM1 in NGC 2264-D, M DM mi ~ l.lxlO -4 M^yr" 1 
(cf § 7). The mass accretion rate onto C-MM3 is thus likely 
in the range l.lxlO -4 < McMMa < 3xl0~ 3 M .yr -1 , 
which is sufficient to overcome the radiation pressure of a 
~ 10 M Q star (Wolfire & Cassinelli 1987, Jijina & Adams 
1996). 

Assuming a typical core star formation efficiency e ~ 50%, 
consistent with calculations of the feedback effect due the 



protostellar outflow (e.g. Matzner & McKee 2000), C- 
MM3, which is thought to be a ~ 40 M Q core, could then 
be the progenitor to a ~ 20 M@ star. 

In the context of the turbulent core model, the for- 
mation of a 20 Mq star requires a velocity dispersion at 
the surface of the initial prestellar core a s ~ 1.0 km.s -1 , 
assuming the fiducial values of MT03 and a clump sur- 
face density S c ; ~ 0.6 g.cm~ 2 . The model then predicts 
a mass accretion rate m* ~ 2.0xl0 -4 M Q .yr _1 when 
half of the final stellar mass has been accreted (i.e., end 
of the Class phase), rising to a maximum value of 
~ 2.8xl0~ 4 M .yr _1 at the end of the accretion phase. 
In this model, the total star formation time from the birth 
of the protostar to the moment it reaches its final mass, 
m+f = 20 M Q , is t*/ ~ 1.7xl0 5 yr. We conclude that the 
observed properties of C-MM3 in NGC 2264-C are roughly 
consistent with, albeit somewhat more extreme/dynamic 
than, the MT03 model. 

The dynamical time estimated in § 6, corresponding 
to the time for the two Class protostars C-MM2 and C- 
MM4 to reach the position of C-MM3, is ~ 1.7 x 10 5 yr. 
This timescale is similar to the formation time of a 20 M© 
star in the MT03 model and to the duration of the 
Class I protostellar accretion phase for low-mass stars (e.g. 
Greene et al. 1994). Therefore, significant evolution evo- 
lution of C-MM2 and C-MM4 is likely to occur before 
a merger can take place in the center of NGC 2264-C. 
Nevertheless, a strong dynamical interaction between C- 
MM2 and C-MM4 seems possible before the associated 
protostars have achieved their final masses. Accordingly, 
we suggest that the formation of an ultra-dense protostel- 
lar core through the merger of two or more intermediate- 
mass protostars in the central part of a protocluster may 
be a plausible route to high-mass star formation. Although 
in itself such a merger does not directly overcome the prob- 
lems arising from radiation pressure, it provides a dynami- 
cal way to build up a massive protostellar core with phys- 
ical characteristics similar to the McKee & Tan model. 
Albeit reminiscent of the Bonnell et al. (1998, 2002) sce- 
nario, this suggested route differs in that it would be based 
on the collision and coalescence of protostellar dense cores 
in a gas-dominated environment, as opposed to stellar 
mergers in a stellar-dominated potential. 

8.3. Conclusions: Massive star formation in clusters 

The picture of massive, clustered star formation which 
emerges from the present study is intermediate between 
the highly dynamic scenario of Bonnell et al. (1998) and 
the turbulent core scenario of MT03. In agreement with 
the first scenario and at variance with the second scenario, 
our observations support the view that protoclusters are 
in a state of global, dynamical collapse rather than ap- 
proximate hydrostatic equilibrium. On the other hand, in 
contrast to the first scenario and better agreement with 
the second scenario, most stars in a cluster, including mas- 
sive stars, appear to acquire their final masses while still 
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embedded in a gas-dominated environment, with processes 
occurring later on playing relatively little role. 
Based on our observations, we propose a mixed scenario of 
protocluster formation according to which a pre-existing 
massive, turbulent clump is strongly compressed and in- 
duced to collapse from outside. Such a clump rapidly en- 
ters a phase of global, near free-fall collapse leading to 
an averaged density profile approaching p cx r -15 in the 
central regions. (Note that the presence of large-scale, su- 
personic inward motions may well be a generic feature 
of all embedded protoclusters in high-mass star-forming 
complexes - Motte et al. 2005.) As the clump is turbulent 
and contains many Jeans masses at the onset of dynami- 
cal collapse, it quickly fragments and produces a number 
of relatively massive 'cores'. These cores which are denser 
than the clump material surrounding them are collapsing 
on significantly shorter timescales than the entire proto- 
cluster clump. Because of the high turbulent /dynamic am- 
bient pressure, the cores are characterized by large densi- 
ties, reminiscent of the MT03 model, and form individual 
protostars with large mass infall rates. Although most of 
these cores have masses exceeding the local Jeans mass, 
we speculate that further fragmentation will be largely 
inhibited as the system is already dominated by strongly 
converging collapse motions as opposed to random tur- 
bulent motions 2 . This contrasts with the collapse of mas- 
sive turbulent cores initially in hydrostatic equilibrium, 
which yields multiple low-mass fragments according to re- 
cent numerical simulations (Dobbs et al. 2005). In other 
words, we speculate that the state of global, dynamical 
collapse of the parent protocluster clump will strongly 
limit the process of sub-fragmentation in the cores, which 
may help to solve one of the main problems of the MT03 
model of massive star formation. In our radiative trans- 
fer model of NGC 2264-C, the 3D turbulent velocity dis- 
persion (7i Mr £>-3D=0.81 km.s -1 is smaller than the infall 
velocity V m/ = 1.3 km.s -1 at r < 3.5 x 10 4 AU, so that 
coherent infall motions dominate over turbulent motions 
at small radii. While the estimated level of turbulence is 
probably still sufficient to alter the gas motions from a 
pure radial infall pattern, we thus do not expect turbu- 
lent, Jeans- type fragmentation to dominate in the inner 
^, 3.5 x 10 4 AU region. On the contrary, our observational 
results rather suggest that two or more cores are in the 
process of merging in this central region. 

We conclude that the millimeter continuum sources 
identified in the central part of NGC 2264-C are the 
probable precursors to intermediate-/high-mass stars. 
Powerful winds will then be likely generated that can clear 
away most of the protocluster gas and thus effectively 



2 Strictly speaking, the observations discussed here mostly 
constrain the motions of the parent clump on scales larger 
than ~ 10 4 AU. Higher-resolution observations with the IRAM 
Plateau de Bure interferometer have been undertaken to probe 
the dynamics and degree of fragmentation of individual MMSs 
or 'cores' on smaller scales and will be reported in a future 
paper. 



stop the global collapse of the system. In this way, the 
protocluster clump NGC 2264-C will likely evolve into a 
revealed star cluster similar to the neighbouring near-IR 
cluster discussed by Lada & Lada (2003) (see Fig. I16|) . 



9. Summary 

The main results of our comprehensive millimeter study of 
the cluster-forming clumps NGC 2264-C and NGC 2264-D 
are as follows: 

1. From our dust continuum maps, we have extracted 
a total of 27 compact millimeter sources (MMSs) in 
NGC 2264, 12 in clump C and 15 in clump D. These 
MMSs have a typical diameter ~ 0.04 pc and masses 
ranging from ~ 2 to 41 M Q . The median mass of the 
MMSs is slightly larger in NGC 2264-C (~ 10 M ver- 
sus ~ 6 M Q in NGC 2264-D), which is the more mas- 
sive clump with a total gas mass ~ 1650 M & (com- 
pared to - 1310 M e for NGC 2264-D). In both cases, 
the mass contained in the MMSs amounts to less than 
10% of the total mass of the parent clump. 

2. The MMSs of NGC 2264 exhibit broad N 2 H+(l-0) 
linewidths tracing supersonic line-of-sight velocity dis- 
persions dime ~ 0.7 km.s -1 ) on ~ 0.1 pc scales. This 
is ~ twice and ~ five times as large as the line-of-sight 
velocity dispersions measured on similar scales toward 
the dense cores of the p-Ophiuchi main cloud and the 
isolated prestellar cores of Taurus, respectively, sug- 
gesting the MMSs of NGC 2264 are significantly more 
turbulent than the dense cores of both p-Oph and 
Taurus. 

3. As much as 70 % of the MMSs of NGC 2264-C are not 
starless but already contain candidate Class proto- 
stellar objects traced by the presence of shocked H2 
jets. In NGC 2264-D, we have direct evidence of the 
protostellar nature of only 25% of the MMSs from 
the close association with 2MASS and/or MSX near- 
/mid-IR source(s). This apparent difference between 
NGC 2264-C and NGC 2264-D may be an artifact 
due to the lack of systematic jet/outflow surveys in 
NGC 2264-D. The high percentage of candidate Class 
objects in NGC 2264-C may result from a short, effi- 
cient episode of star formation, as expected in the case 
of triggered cloud collapse. 

4. In NGC 2264-C, we observe widespread blue-skewed 
"infall" line profiles in optically thick tracers such as 
HCO + (3-2) or CS(3-2), consistent with the presence 
of large-scale infall motions in this clump as previously 
suggested by Williams & Garland (2002). 

5. Furthermore, our mapping in low optical depth tracers 
such as N 2 H + (l-0) or H 13 CO + (1-0) reveals a new ve- 
locity feature in the central part of NGC 2264-C, in the 
form of a sharp velocity discontinuity ~ 2 km.s -1 in 
amplitude centered on the millimeter source C-MM3. 

6. Taken together, we interpret these two velocity fea- 
tures as the signatures of the large-scale collapse of a 
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cylinder-like, filamentary structure along its long axis. 
Radiative transfer calculations confirm this view and 
suggest an infall velocity V in f—1.3 km.s -1 along the 
main axis of the NGC 2264-C prolate clump, corre- 
sponding to a total inflow rate ~ 3x 1CT 3 M Q .yr _1 
toward the central Class object C-MM3. 

7. We have also found evidence for both infall and rota- 
tion motions toward the Class object D-MM1 in the 
NGC 2264-D clump. Based on radiative transfer mod- 
elling of our HCO + and CS observations, we estimate 
an infalling velocity Vj n /=0.1 km.s -1 which, given the 
high density of D-MM1, corresponds to a mass infall 
rate ~ 1.1 x 10~ 4 M Q .yr _1 . Such an infall rate is in 
principle high enough to overcome the radiation pres- 
sure of a massive }t 10 M star. 

8. Comparison of our observations with existing scenar- 
ios for massive star formation suggests that we are 
witnessing the initial stages of the formation of a high- 
mass ~ 10-20 M© star in the central part of NGC 2264- 
C. Based on this comparison, we propose a picture of 
high-mass star formation in protoclusters intermedi- 
ate between the highly dynamic scenario of Bonnell et 
al. (1998) and the turbulent core scenario of McKee 
& Tan (2003). In this observationally-driven picture, 
the large-scale, dynamical collapse of a massive, ini- 
tially unstable clump would lead to the formation of a 
turbulent, massive, and ultra-dense core with proper- 
ties reminiscent of the McKee & Tan model, through 
the gravitational merger of two or more intermediate- 
mass Class protostellar cores in the central part of 
the clump, while the potential of the system is still 
gas-dominated. 
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Fig. 6. N 2 H+(101-012) and H 13 CO+(1-0) position-velocity diagrams observed in the central part of NGC 2264-C along 
four axes shown schematically on the left. From top to bottom: the first row corresponds to PV diagrams taken along 
an East- West (EW) axis going through C-MM2, C-MM3, and C-MM4; the second row corresponds to a North-South 
axis going through C-MM3 (NS-3); the third row is for a North-South axis going through C-MM4 (NS-4); the fourth 
row is for a North-South axis going through C-MM2 (NS-2). Contour levels are 1.1 and 1.8 K for the N 2 H+(101-012) 
diagrams, 1.05 and 1.4 K for the H 13 CO + (1-0) diagrams. 
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Fig. 10. HCO+(3-2), H 13 CO+(3-2), H 13 CO+(1-0) (top) and CS(3-2), CS(5-4), C 34 S(2-1), C 34 S(3-2), C 34 S(5-4) (bot- 
tom) spectra observed along an East- West axis going through C-MM2, C-MM3, C-MM4 in NGC 2264-C. Synthetic 
spectra corresponding to the "best-fit" 2D radiative transfer model described in the text (see Fig. IHt for input param- 
eters) are superimposed in light grey. The central position corresponds to C-MM3. 



